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We propose how to create a non-trivial Haldane phase in atomic two-component Fermi-gas loaded 
on one-dimensional (1-D) optical lattice with trap potential. The Haldane phase is naturally formed 
on p-band Mott core in a wide range of the strong on-site repulsive interaction. The present proposal 
is composed of two steps, one of which is theoretical derivation of an effective 1-D S = 1 interacting- 
chain model from the original tight-binding Hamiltonian handling the two p-orbitals, and the other 
of which is numerical demonstration employing the density-matrix renormalization-group for the 
formation of the Haldane phase on p-band Mott core and its associated features in the original 
tight-binding model with the harmonic trap potential. 
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Ultra-cold Fermi gas loaded on optical lattice (FGOL) 
is one of the most fascinating systems in modern physics 
[l|. Its controllability is far beyond our past image and 
experience in condensed matters. For instance, the two- 
body interaction between Fermi atoms is widely variable 
from strongly attractive to repulsive. Moreover, filling 
and population imbalance are freely tunable, and lat- 
tice geometry is highly flexible. Thus, FGOL is regarded 
as an outstanding simulator for strongly-correlated elec- 
tronic materials, in which a great number of controversial 
issues like high-temperature superconducting mechanism 
in cuprate superconductors still remain unsolved. 

Very recently, there have been intense theoretical and 
experimental interests in extending the energy band pro- 
duced by the optical lattice (OL) potential from the 
ground single one to higher multiple-ones by utilizing 
multi-degenerate higher orbitals Such multiple or- 

bital degeneracy allows us to study orbital degrees of free- 
dom in addition to charge and spin ones. Then, FGOL 
becomes a more realistic and rich simulator, in which 
orbital ordering and high-spin correlation are accessible 
subjects like real solid-state matters including transition 
metals and other heavy elements 0, H|. In this paper, 
we study one of the most fundamental issues associated 
with the multiple band degeneracy. The present target 
is formation of a gapped quantum phase using double p- 
orbital degeneracy under a one-dimensional(l-D) system, 
as shown in Fig.l. We show that an effective low-energy 
Hamiltonian is given by the S = 1 Heisenberg model and 
the gapped Haldane phase [6] is widely sustained by the 
harmonic trap potential to confine Fermi atoms. The fer- 
tility of the Haldane phase, i.e., gapful spin-excitation, 
nonlocal string order, and spin- 1/2 edge magnetization 
induction |7l-[lQ| can be systematically explored because 
of the wide controllability in FGOL's. 

Feasibility studies to create the Haldane phases in 
atomic gas OL's have been made on spin 1 bosonic gases 
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FIG. 1: (Color Online) Schematic figures for a 1-D Fermi 
gas. (a) An optical lattice (OL) potential along ^-direction 
and a cylindrically-symmetric vertical trap-potential inside 
x;?/-plane. The latter causes to form discrete levels as de- 
picted at the right panel, in which the lowest level is s-orbital 
while higher two degenerate ones are p^y-obitals. The Hal- 
dane phase is formed under fully-filled s- and partially-filled 
p-orbitals (see text), (b) p^-orbitals formed inside the ver- 
tical trap and their p-bands creatd by the OL. When p-bands 
are active, this system is described by a multi-band Hub- 
bard chain with the on-site intra-orbital interaction (U pp ) and 
inter-orbital one (U PxPy ) (see Eq.© and text). 



loaded on 1-D OL with one atom per site [TTJ] , two- 
component 1-D FGOL utilizing 1st and 2nd Bloch-bands 
with Feshbach resonance [12], and ultra-cold fermionic 
171 Yb in 1-D OL . The Haldane insulator, in which 
a charge gap opens according to Haldane's conjecture 
was also predicted for bosonic gases with dipole inter- 
action [Ullll and multi-component 1-D FGOL [ItJ. In 
this paper, through the derivation of an effective 5 = 1 
interacting-chain model and the numerical simulations of 
a multi-band Hubbard model, we show that the Haldane 
phase is naturally formed in a two-component 1-D FGOL 
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with the emergence of p-band Mott core. We also reveal 
features of the Haldane phase in trapped systems. 

We start with the following Hamiltonian describing a 
1-D elongated FGOL schematically depicted in Fig.l, 



H 



E 



dx 



(i) 



with ho = (—h 2 /2m)V 2 + V ver + V op t and the coupling 
constant of the two-body interaction g. The cylindrically- 
symmetric vertical trap (on xy-plane) and the OL poten- 
tial (along z-axis) are, respectively, V ver and V op f 

A multi-band Hubbard-type model is derived from 
Eq.(l) using the expansion ip a = ^2 a ^ J2i c a,(3,<j,iUaWp,i , 
where u a and wpj are a wavefunction associated with 
the eigensystem of [(— h 2 /2m)\/ 2 L + V ver ] u a = e a u a and 
a Wannier function formed by the OL potential. The in- 
dices a and f3 represent discrete levels caused by the trap 
potential Ker and Bloch-band by the OL potential, re- 
spectively. Here, Ker is not so tight that the 2nd p-orbital 
energy-level does not exceed over that of 2nd Bloch-band 
and the 2nd level is partially filled, as shown in Fig. 1(b). 
Hereafter, ignoring the contribution from higher Bloch- 
bands, we drop the index f3 . Then, including the 2nd 
level corresponding to p x ( y ) -orbital and taking the tight- 
binding approximation, we obtain a multi-band 1-D Hub- 
bard Hamiltonian [l8| 



H ~ ^1 ^1 t c lx,(T,i c oc,crJ + ^1 e a n a,crA 
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where n a ,a y i(= (ji c OL ^^i) is the on-site number-density 
operator of a-orbital with pseudo-spin a and S a ,i = 
\ a' c a i C r r cr,cr / c Q! ,i,cr / with Pauli matrices r is the lo- 
cal spin- 1/2 operator. The summation of a is taken 
from the ground s-orbital to the 2nd degenerate p- 
orbitals (p x and p y ) as shown in Fig. 1(c). t and U aa ' 
are the hopping and on-site interaction energy inte- 
grals defined as t = - J dzw i+1 [^^2 + K P t) Wi and 

Uaa' — 9 J dxw\v? a VL 2 ot , , respectively. The in-plane cylin- 
drical symmetry of u a (x±) gives the relations U sPx — 
U sPy (= U sp ) and U PxPx = U PyPy (= U pp ). A significant 
feature in Eq.(j2j) is the existence of Hund's like terms 
— Uaa' Sa a ' So,' a , which captures the so-called orbital 
physics. Similar Hamiltonian to Eq.(j2j) may also ap- 
pear in ultra-cold fermionic 171 Yb with 1 5'o- and 3 Po- 
states [H [l9(| . 

Now, we derive aS = l Heisenberg chain from Eq.(j2j) 
in the strong coupling limit (U aa / ^> t). We con- 
sider a case when the 1st orbital is fully occupied and 
2nd ones are half-filled. Now, this system seem to be 



similar to two Heisenberg chains coupled by ferromag 
netic exchange interaction, which can produce the Hal- 
dane phase I20l . In fact, the second order perturba- 
tion scheme [21j leads to Hj = J ex ^ 



= 2t 2 /(U pp 
= ±Y 

2 ^-^a—p x , 



^PxPy 



Ki j> Si - Sj, where 
) and the local spin-1 operator 



rC a ,a,i [18]. Hence, 



s l 

gapped Haldane phase can emerge for large U pp if the 
2nd p-levels are half-filled. In FGOL experiments, this 
requirement is achievable in a central Mott core formed 
by interplay between a trap potential and large U pp . 

Next, let us turn to the numerical demonstration of the 
Haldane phase by the density-matrix renormalization- 
group (DMRG) method [H, [23|. Typical atomic-gas ex- 
periments employ harmonic trap potential in all direc- 
tions to prevent escape of atoms. We add the harmonic 
trap along ^-axis, V ho (i) = V[2/(L - l)] 2 [i - (L + l)/2] 2 
to our system, as well as the aforementioned vertical trap 
V^er- Here L means the total number of sites. The total 
Hamiltonian is H + a i Vh (i)ria,crA' Again, we drop 
off the terms in the Hamiltonian relevant with the 1st 
orbital fl8l|. We simulate the model handling only 2nd 
p-orbitals [24]. In the 2-band Hubbard model, the on-site 
intra-orbital interaction U pp and inter-orbital interaction 



U PxPy are mutually connected via U Pa 



9 u pp- 



We 



note that U aa ' = 9 J dxwfu^u 2 ^, [25|. 

Now, we present numerical results. We obtain 
particle-density distribution n(i) = J2 a = 



=Px,Py, 



-P X ,Py 



(n, 



OLy\ A ' 



and a spin-density distribution m(i) = J2 a= 
n a,l,i) using DMRG simulations. Figure [2fa) shows n(i) 
for Upp/t = 15 and V/t = 10. Owing to Vh o (0? P~ 
band Mott plateau is formed in the trap center and sur- 
rounded by regions with trap-gradient-dependent filling 
(below half). Varying an input imbalance ratio p = 
rrii/ni, which is initially set in experiments, we mea- 
sure the normalized polarization on the Mott core [26| , 
M = X^eMott m i/ n ii as seen in Fig.[2fb). For compari- 
son, we also show in the inset the result for U PxPy = 0, 
which corresponds to two identical 1-D Hubbard chains 
without coupling. A plateau-like behavior occurs on the 
Mott-core polarization up to a critical imbalance ratio 
p c (~ 0.17) as expected in the presence of the gap, while 
for U PxPy = the Mott-core is smoothly magnetized. A 
gapful phase for the spin flip excitation is found in the 
p-band Mott-core region. Thus, one of the key features 
of the Haldane phase is shown. 

Let us more carefully examine the trap center and 
the outer regions for a balanced case p = (5^(= 
^ZiS^) = 0) and a slightly-polarized case p = 0.0125 

(Stot = 1). First, we focus on the spin-density dis- 
tribution profiles on the p-band Mott core. Figure [3] 
shows that staggered magnetization structures emerge 
at the edges of the p-band Mott core and exponen- 
tially decay toward the trap center. This structure is 
the same as the typical one observed at open boundary 
edges in S = 1 Heisenberg chains 0, [Toj| and is known 
to be induced by S = 1/2 fraction left at open bound- 
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FIG. 2: (Color Online) (a) A typical spatial distribution pro- 
file of the particle density n(i) = J2&= Px Py So- n ^,a,i for bal- 
anced 160 fermions (80 t? 80 \) with the interaction constants 
Upp/t — 15 , Up xPy — \Upp , and the trap potential strength 
V/t — 10. (b) The Mott-core polarization M vs. the popu- 
lation imbalance ratio p. The employed parameters are the 
same parameter as in (a) except for p. For comparison, the 
inset is the case of U pp /t = 15 and U PxPy /t = 0, which is 
equivalent with two independent S = 1/2 interacting chains. 
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ary edges. In the present case, the p-band Mott core 
is not sharply terminated by the open boundary con- 
dition. The edge magnetization amplitude in the in- 
duced staggered structures is sizeably reduced, compared 
to the results of S = 1 Heisenberg chain in the open 
boundary case [9] . This difference comes from coupling 
between S = 1/2 fraction left at the edges of the p- 
band Mott core and a ferromagnetic metallic phase in 
the outer regions. The ferromagnetic property in the 
outer regions is confirmed by calculating spin-gap energy 
AE P = E p — E p =o, polarization on outer regions, and 
the Mott-core polarization [l8j . Here, E p is the ground- 
state energy for the population imbalance ratio p. This 
ferromagnetism is also consistent with theoretical predic- 
tion for a uniform hole-doped two-degenerate-band sys- 
tem with strong Hund's coupling [27j. Next, let us exam- 
ine the outer regions to understand coupling between the 
ferromagnetic metal and the 5 = 1/2 fraction at the edge 

of the Mott core. As for S^l = 0, as shown in Fig.[3^a), 
the spin-density profile is anti-centrosymmetric. We cal- 
culate the integral of m(i) over the left (right) outer re- 
gion, m L(R) EE EiGleft (right) 771 W' Then > We find that the 

left and right regions, respectively, are polarized as one 
up-spin (mL ~ 1) and one down-spin (mR ~ —1). For 

Sfol = 1, as shown in Figj3^b), both regions are polarized 
as one up-spin (mL = mR ~ 1). We always find similar 
features for various p. The polarization difference be- 
comes either |mL — raR,| — or |mL — ttir\ ^ 2 as shown in 
Fig-El^c). These results mean that S = 1/2 fractions left 
at the edges of the p-band Mott core may spread over the 
outer regions. We further characterize the ferromagnetic 
metal on the outer regions. We compare the spin-gap 
energy with Mott-core polarization. The ground-state 
degeneracy (i.e., AE P = 0) exists up to p ~ 0.11, as seen 
in Fig.[3jd), while plateau-like behavior in Fig.[2fb) still 
continues up to p ~ 0.17. This result implies that all 



FIG. 3: (Color Online) Distribution profiles of particle den- 
sity (dashed blue line) n(i) and spin density (solid red line) 

m W = J2 a =p x , Py K,t,i-Wa,i,i) for spin imbalance (a) p = , 

(Sfcl = 0) and (b) p = 0.0125 , (Sfcl = 1) . The inset of (a) 
is a focus on the spin density around the edge of the left 
Mott core, (c) The polarization in the left (right) surround- 
ing phase ra L (R) = Z)ieieft(right) m W versus the population 
imbalance ratio p. (d) The spin gap AE P = E p — E p =o versus 
the population imbalance ratio p . All other physical param- 
eters are the same as in Fig|2ja). 



spin degrees of freedom in the system accumulate on the 
ferromagnetic metal before p ~ 0.11 . Therefore, the fer- 
romagnetic metal on the outer regions corresponds to an 
edge state in uniform S = 1 Heisenberg model. 

Here, we mention the case of small Upp/t. The Haldane 
phase may survive in small U region as discussed in the 
case of multi-component 1-D FGOL [ItJ • We emphasize 
that, for a wide-range of the interaction strength and the 
trap potential strength V/t, the gapped Mott core more 
or less appears and the Haldane phase is realized in the 
central core region [l8[ . 

We further increase the population imbalance ratio to 
examine the spin-density distribution profiles before and 
after p c . As for p = 0.1625 (N^ — N± = 26) , as seen in 
FigHfa) , the central Mott core is not magnetized, and all 
the polarization concentrates on the ferromagnetic metal- 
lic phases, while the magnetization amplitude is consid- 
erably enhanced compared to Fig|3l A coupling with the 
Haldane phase and polarized ferromagnetic metal brings 
about the enhancement of the edge magnetization. On 
the other hand, when p = 0.475(> p c ), the staggered 
structure disappears, and even the p-band Mott core is 
uniformly magnetized, as shown in FigH^b). As for an 
extremely population imbalance ratio p = 0.9375, a pe- 
riodical oscillation appears inside the p-band Mott core 
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FIG. 4: (Color Online) Particle density (dashed blue line) and 
spin density (solid red line) distribution profiles for population 
imbalance ratios (a)p = 0.1625 (93 f, 67 I), (b)p = 0.475 
(118 t,42 I) , and (c) p = 0.9375 (155 t, 5 I). All other 
parameters are the same as in Fig[2] (a). 

as shown in FigH^c). This oscillation can be explained 
as follows. A low-energy excitation of the spin- S Heisen- 
berg model under strong magnetic field is describe d by 
the Luttinger liquid irrespective of the spin length S [28[ . 
Then, this theoretical picture predicts that the period- 
icity of a spin-density wave (SDW) is characterized by 
2fcp, where fcp is the Fermi wave vector in the equiv- 
alent spinless fermion system. In the present trapped 
systems, the oscillation periodicity is g iven by the rela- 
tion 2k F = (27r/L Mo tt)A r Mott,; 0, Ha, where L Mo tt is 
the size of the Mott core, and iVMott,4, * s tne total num- 
ber of spin-down component inside the Mott-core. In- 
deed, the magnetization profile shows an expected oscil- 



lation. In Fig.|4^c), the number of the minimum peaks 
of the SDW is 5. This number should be almost equal to 
^Mott,|- In fact, we can directly evaluate iVNott,^ and 
find A^Mott,! — 5. A more systematic analysis on p- 
dependence in such a SDW-like oscillation is shown in 
Ref. [H|. 

In conclusion, we confirmed that the Haldane phase 
emerges as a leading phase when utilizing 2nd degener- 
ate p-orbitals in 1-D FGOL. At the half- filling condition 
for p-orbitals in the large U limit, the multi-band Hub- 
bard Hamiltonian was reduced to the 5 = 1 Heisenberg 
chain model. Then, the emergence of the Haldane phase 
and its associated physics were demonstrated by DMRG 
studies on the original model together with the harmonic 
trap potential. The polarization on the outer regions is 
easily observed in experiment [50, [H| and spin-selective 
single-site addressing [32] can be used to detect the stag- 
gered magnetization on Mott core. The formation of the 
Haldane phase in 1-D FGOL allows us to not only inves- 
tigate its features but also open a new avenue towards 
topologically-protected quantum state engineering [33[ • 
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